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模糊 数学 评价 河流 水 质 的 模型 探讨 
孙 划 
(湖南 省 环境 监测 中 心 站 ， 长 沙 410004) 


摘 要 采用 模糊 数学 中 模糊 聚 类 的 2 种 方法 即 最 大 给 阵 元 法 和 模糊 相关 分 
析 法 对 武 水 河水 质 进行 分 析 、 聚 类 和 评价 ,试图 通过 实例 阐明 最 大 算 阵 元 法 和 
模糊 相关 分 析 法 对 河流 水 体 单 元 、 尤 其 是 污染 参数 进行 聚 类 分 析 评价 ,从 不 同 
侧面 全 面 剖 析 河 流水 质 污染 状况 ， 分 析 污 染 因子 在 河水 中 的 相互 作用 及 各 种 
物理 化 学 作用 所 形成 的 污染 参数 之 间 的 相互 关系 。 

关键 词 水 质 模糊 聚 类 分 析 评价 


水 质 评价 是 环境 科学 中 广泛 使 用 的 一 种 
表述 水 体 水 质 状况 的 手段 ， 其 评价 方法 目前 
通用 的 是 单项 污染 指数 法 。 单 项 指数 法 作为 
一 种 标准 方法 具有 简明 操作 性 强 等 特点 ,但 
为 了 更 深入 全 面 地 分 析 河 流水 质 ， 不 妨 在 水 
质 评价 方法 学 上 进行 更 多 的 探讨 。 

水 体 污染 程度 和 污染 因子 在 河流 中 的 相 
互 关系 等 都 客观 存在 着 模糊 性 。 用 模糊 数学 
中 模糊 聚 类 的 2 种 方法 从 不 同 角度 对 水 体 水 
质 进 行 痢 析 ， 对 水 体 单 元 和 污染 因子 进行 聚 
类 分 析 ， 试 图 通过 这 2 个 模型 ， 建 立 用 模糊 
聚 类 分 析 对 水 体 单 元 和 污染 物 进 行 聚 类 评价 
的 通用 方法 。 


1 模糊 数学 进行 水 质 评价 的 理论 和 模型 


1.1 模糊 矩阵 及 其 建立 

定义 矩阵 R= (rw 叫做 一 个 模糊 矩 
阵 、 对 任意 i =1,2,…n, 7 二 1, 2,，…m， 者 
有 rE 区 0 

用 模糊 数学 进行 水 质 评价 ， 建 立 模糊 和 拖 
阵 是 一 个 核心 工作 ,对 于 任何 一 个 研究 水 体 ， 
都 有 两 个 因素 构成 ， 即 水 体 单 元 和 各 水 体 单 
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元 中 的 污染 参数 。 其 模糊 相关 和 矩阵 的 矩阵 元 
”的 确定 如 下 : 

设 X= (rzr，…z) 为 待 分 类 事物 的 
全 体 ， 即 水 体 单元 ，zi 由 一 组 数据 U,,U,;， 
…U7w 表 征 。 和 上 的 模糊 相关 答 阵 民 一 (mn 
可 用 多 种 方法 计算 。 在 此 ， 为 建立 两 个 水 体 
单元 之 间 的 相似 关系 ， 采 用 模糊 集中 的 “ 贴 
近 度 ”概念 之 一 的 最 大 最 小 法 计算 : 

5 min (Ure 大 

一 二- (1, Jn) (1) 


7,) 


max (Ui;,, LU) 
上 式 中 水 质 因子 对 目标 值 隶 属 度 Lo 的 








确定 按 下 列 公 式 计算 : 
we 这 Cu>S， 
1 ”实测 值 优 于 目标 值 时 
式 中 : U, 一 一 i 水 体 单元 ; 类 污染 物 对 目标 
值 的 隶属 程度 ; 
C, 一 一 i 水 体 单元 j 类 污染 物 的 实测 
浓度 ; 
S, 一 一 i 水 体 单元 j 类 污染 物 的 目标 
值 。 


1.2 污染 因子 的 模糊 相关 分 析 
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对 于 环境 问题 ， 在 引起 环境 变化 的 诸 因 
素 中 ,往往 有 未 知 的 或 无 法 监测 的 ， 就 是 在 
已 监测 的 因子 中 ， 很 可 能 有 些 相 互 没有 实质 
性 的 联系 或 由 于 某 种 干扰 使 得 有 实质 性 联系 
的 因素 未 显示 出 来 。 因 此 有 必要 采取 一 定 的 
数学 方法 进行 筛选 。 

以 相似 系数 法 建立 模糊 相关 矩阵。 由 此 
矩阵 可 分 析 各 样本 间 的 相互 关系 ， 即 各 污染 
因子 间 的 相关 性 。 但 应 注意 的 是 ， 此 处 的 相 
关 是 两 个 污染 因子 之 间 在 量 上 的 相互 依赖 程 
度 ， 至 于 这 种 依赖 是 正 相 关 还 是 负 相 关 则 需 
要 依 其 物理 化 学 机 制 而 定 。 

在 已 建 的 目标 值 隶 属 度 的 基础 上 ， 以 相 
似 系 数 法 建立 模糊 相关 矩阵 ， 其 矩阵 元 x， 


0. 96 1 

0. 98 0. 94 ] 
R= 10.99 0. 96 0. 98 

0. 98 0,. 95 0. 98 

0. 94 0. 91 0. 93 

Qs95 O05:93 0. 95 


以 上 矩阵 为 模糊 相似 矩阵 ， 模 糊 相 似 敌 
阵 具有 自 反 性 和 对 称 性 ,但 不 具有 传递 性 ,只 
有 将 模糊 相似 关系 通过 求 其 传递 闭 包 而 转换 
为 模糊 等 价 关系 , 扰 阵 的 任 一 截 和 矩阵 Ri 才 具 


0. 96 1 

0. 98 0. 96 1 
R'= ROR’:= |0. 99 0. 96 0. 98 

0. 98 0. 96 0. 98 

0. 94 0. 94 0. 94 

0. 95 0. 95 0. 95 


可 将 其 中 的 元 素 归 并 分 类 ， 当 和 从 1 降 为 0， 
由 于 R; 不 断 变化 ， 分 类 由 细 变 粗 逐 渐 归 并 ， 


0. 94 
0, 95 


为 : 
pO 
k=1 


> (Us—U,)? | (Ui—U,) 
二 1 太一 1 


= ] a 加 sn 
式 中 . LU 一方 Zu U ,= 一 ZU 


(2) 


Ti; 





2 应 用 实例 


2.1 最 大 和 矩阵 元 法 对 武 水 河水 质 的 评价 
2.1.1 水 体 单 元 的 聚 类 评价 

以 监测 断面 为 水 体 单元 ， 根 据 各 水 体 单 
元 水 质 监测 值 建立 各 水 体 单元 的 水 质 因 子 对 
环境 目标 的 隶属 度 ， 由 隶属 度 求 得 的 模糊 相 
关 和 矩阵 为 : 


1 (3) 
0.98 | 
0. 94 
0. 95 0. 93 ] 


有 传递 性 ， 即 有 尺 为 和 一 截 矩 阵 《〈 布 尔 拖 
阵 ) ,经 过 R 和 矩阵 的 变换 ,得 RR 的 等 价 和 矩阵 如 
下 : 

其 中 4 为 置信 水 平 。 对 于 一 个 等 价 和 矩阵 


1 (4) 
0. 98 lL 
0. 94 0. 94 1 
0. 95 0. 95 0. 94 1 
形成 动态 分 类 。 


当 取 不 同 的 置信 水 平时 ， 各 水 体 单元 得 


。25 。 


2000 年 第 4 期 贵州 


到 分 类 。 其 依据 是 各 水 体 单元 受 污 染 程度 的 
相似 性 及 差异 性 。 当 置信 水 平 越 高 ， 水 体 单 
元 得 到 越 详细 的 分 类 ; 当 置 信 水 平 越 低 ， 水 
体 单元 分 类 的 模糊 性 越 大 。 当 到 和 = 0.98、 
0. 96 时 武 水 河 宜 章 江 段 7 个 水 体 单元 的 分 
类 如 下 ， 
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[at | 
[LA]; 

[A As], [As、 A Aj, LAs], Lj, 

武 水 河 宜 章 江 段 水 体 单元 动态 分 类 聚 类 
如 图 1 所 示 。 同 时 ， 可 根据 相似 矩阵 中 每 行 
最 大 值 即 置信 水 平 , 将 水 体 单 元 分 级 ( 表 1) 


[4 ， LA;.、 A1、4A4sj]， Ee 


\ ls 3 A; A /7 分 类 数 
(0.99 
0.98 5 
0.96 4 
0.95 ) 
(}.94 I 

图 1 武 水 河 宜 章 江 段 水 体 单元 动态 分 类 聚 类 
表 1 武 水 河 宜 章 江 段 水 体 单元 的 分 级 
序号 1l 2 2 + 5 6 7 
水 体 单 元 和 Ai A Als 1 4 /ly As 

入 0. 99 0. 99 们 .98 0. 90 (站 ,96 0.95 0. 9] 

分 类 1 


根据 以 上 分 析 看 出 ， 武 水 河 宜 草 江 段 7 
个 断面 其 污染 程度 的 排序 为 4; 二 4A; 二 4, 
4，45>4、 4 第 1 第 4 两 个 断面 污染 程 
度 相 近 ， 属 于 受 污 染 最 轻 的 断面 ; 第 6 断面 
位 于 两 个 氮肥 厂 排 污 口 下 游 300 m 左右 ， 因 
此 是 污染 最 严重 的 一 个 断面 ; 第 7 断面 虽 距 
第 6 断面 有 较 长 的 距离 ， 得 到 一 定 自净 ， 但 
由 于 期 间 又 有 偏 溪 河 这 一 较 大 排污 支流 的 汇 
和 ,使 这 种 自净 作用 被 抵消 一 部 分 .因此 . 其 
污染 程度 仅 次 于 第 6 断面 ; 污染 程度 位 居 第 
3 的 为 第 2 断面 ， 这 是 因为 第 2 断面 上 游 
200 m 即 为 杨 家 河 排污 支流 ， 杨 家 河 是 水 量 
较 大 ， 污染 较 重 的 排污 口 ， 因 此 使 第 2 断面 
污染 相应 加 重 。 第 2 断面 距 第 4 断面 有 一 段 
较 长 的 距离 。 其 间 虽 也 有 排污 ， 由 于 其 水 量 
和 污染 物 浓度 不 大 .， 因 此 自净 作用 超过 排污 
” 作用， 使 水 质 有 明显 好 转 。 
. 20 四 


采用 式 (1) 所 列 最 大 最 小 法 公式 计算 模 
糊 相 似 关 系 ( 表 2)。 

将 表 2 的 模糊 相似 短 阵 转换 为 模糊 等 价 
和 矩阵 ( 式 (5))。 

当 取 不 同 的 置信 水 平 * 时 ， 各 污染 物 依据 其 
污染 程度 的 大 小 得 到 聚 类 。 

当 A 从 1 取 到 0， 污染 物 的 污染 程度 得 
到 由 细 到 粗 的 分 类 ， 其 分 类 表现 为 一 个 动态 
过 程 ， 各 污染 物 的 聚 类 如 下 : 

4 二 1]， 根 据 污染 程度 污染 物 被 聚 为 7 
类 ， 即 [DO、BOD,、COD、NO, 一 N、NO; 
=Ns CN Pb, Zn)sy [SS CNH= Ns 
[@ 一 OHJj; LAsj; [Fj」; [细菌 总 数 j]。 

4 一 0. 90， 污 染 物 被 聚 为 3 类 ， 即 [DO、 
BOD;,、 COD, NO,—N., B—OH., CN-. As.、 
Pb、Zn、SS、F-]; [NH 一 N]; [细菌 总 
数 ]。 
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表 2 污染 因子 的 模糊 相似 矩阵 
污染 因子 1 2 3 4 5 6 8 9 10 11 12 113 14 
1 DO 1 
2 BOD; ] 1 
3 COD 1 1 1 
| SS 0.93 0.93 0.93 1 
5 NHs—N 0.0!1 0.01 0.04 0.04 1 
6 NO,—N 1 ] 1 Vo O00 
7 NOA—N 1 ] 1 0.93 0.01 1 1 
8 和 -OH 0.94 0.94 0.94 0.89 0.01 0.91 0.94 1 
9 CN- 1 ] 1 9 
10 F- 0.87 0.87 0.87 0.90 0.04 0.87 0.87 0.85 0.87 1 
11 As 0.94 0.94 0.94 0.90 0.04 0.94 0.94 0.90 0.94 0.89 1 
]2 Pb 1 1 1 0.93 0.04 1 国有 
13 Zn 1 T0090 1 0.94 1 0.87 0.91 1 ] 
14 ”细菌 总 数 0.26 0.06 0.26 0.28 0.15 0.26 0.26 0.26 0.27 0.30 0.27 0.26 0.26 1 
| | 1 0.93 0.15 1 ] 0.94 1] 0.90 0.94 1 1 0.30 
| 1 1 0.93 0.15 1 1 0.94 1] 0.90 0.94 ] 1 0.30 
1 1 ] 0.93 0.15 1 1 0.94 1 0.90 0.94 1 1 0.30 
0.93 0.93 0.93 ] 0.15 0.93 0.93 0.93 0.93 0.90 0.93 0.93 0.93 0.30 
0.15 0.15 0.15 0.15 1 HI 0015 O05 .005 O18 0. 1s -00718 015. .0.5 
] 1 ] 0.93 0.15 1 ] 0.94 ] 0.90 0.,94 1 1 .a0 
1 1 1 0.93 0.15 1 1 0.94 1] 0.90 v.94 ] 1 0.30 
人 一 | 4 (5) 
0.94 0.94 0.94 0.93 0.15 0.94 0.94 ] 0.94 0.90 0.91 0.94 0.94 0.30 
1 1 1 0.93 0.15 1 ] 0.94 1 0.90 0.94 ] 1] 0.30 
0.90 0.90 0.90 0.93 0.15 0.90 0.90 0.90 0.90 ] 0.90 0.90 0.90 0.30 
0.94 0.914 0.94 0.93 0.15 0.94 0.94 0.94 0.94 0.90 1 0.94 0.94 0.301 
1 1 ] 0.93 0.15 ] 1 0.94 1 0.90 0.94 1 1 0.30 
1 1 1] 0.93 0.15 1 ] 0.94 1 0.90 0.94 1 1 0.30 
Lo 0.30 0.30 0.30 0.15 0.30 0.30 0.30 0.30 0.30 0.30 0.30 0.30 ] 


根据 以 上 分 析 可 以 看 出 ,污染 物 的 聚 类 ， 
不 仅 反 映 了 污染 物 问 的 相互 依存 关系 ， 即 蘑 
一 污染 物 程度 取决 于 或 依赖 于 其 他 一 种 或 几 
种 污染 因 于 ， 而 且 也 反映 了 哪 一 类 污染 源 对 


水 体 的 作用 最 大 。 


2. 2 污染 参数 的 模糊 相关 分 析 

根据 式 (2) 计算 各 污染 因子 的 模糊 相关 
矩阵 ( 表 3)。 由 表 3 看 出 , DO、BOD;、COD、 
NO.—N.、. NO;—N., CN-. Pb Zn, B®—OH. 


As、F 之 间 相 关系 数 在 0. 98 以 上 , NH 一 N 
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表 3 污染 因子 的 相关 性 


污染 因子 DO BOD, COD SS DN 
D0 

BOD， 和 

COD 1l 1 


SS 0.98 0.98 
Nii=N ‘03 O33. 038 039 1 


NO,—N 1 1 1 0.98 0.33 1 

NO ,一 NM 1 1 1 0,98 0.33 1 
DB--OH 0.98 v.98 0.98 0.96 0.38 0.98 

CN 1 1 1 0.98 0.33 1 
i 0.97 0.97 0.97 0.99 0.43 0.97 


As 0.99 0.99 0.99 0.97 0.38 0.99 


Yo-OHCN- F- As Pb Zn 细 风 


0.98 | 


0.97 0.97 0.97 1 
O99 00 099 "098 l 


Pb 1 1 1 0.98 0.33 1. 0.98 1 0.97 0.99 1 
Zn 1 1 1 0.98 0.33 1 0.98 1 0.97 0.99 1 1 
细菌 总 数 0.52 0.52 0.52 0.58 0.94 0.52 0.52 0.57 0.52 0.64 0.57 0.52 0.52 1 


与 其 它 污染 物 相关 性 较 弱 ， 但 与 细菌 总 数 的 


相关 性 较 强 。 以 上 污染 物 的 相关 性 ， 既 反映 
了 各 类 污染 物 的 来 源 ， 又 反映 了 各 苇 当 因子 
在 水 体 中 的 相互 依赖 程度 。 


3 结果 与 讨论 
3.1 通过 用 最 大 矩阵 元 法 对 水 体 单元 和 汽 


染 因子 进行 聚 类 分 析 ， 同 时 结合 相关 分 析 对 . 


污染 因子 间 的 关系 进行 初步 探讨 ,可 以 看 出 : 
由 于 各 类 污染 物 的 相互 依赖 关系 不 同 及 其 来 
源 不 同 ， 导 致 污染 程度 的 不 同 ， 可 以 通过 模 
糊 聚 类 分 析 这 一 数学 过 程 得 到 体现 。 

3.2 从 对 污染 因子 聚 类 的 结果 看 ,对 水 体 中 
的 污染 物 不 能 孤立 地 看 。 由 于 每 个 污染 因子 
都 具有 一 定 的 来 源 ， 同 时 污染 因子 彼此 间 存 
在 物理 化 学 作用 而 建立 起 的 相互 依赖 关系 ， 
使 得 水 体 中 污染 参数 具有 类 别 的 优势 和 弱 
势 。 就 武 水 而 言 ，NH4 一 N 和 细菌 总 数 占 优 
势 。 
3.3 通过 污染 物 的 聚 类 分 析 ,不仅 揭 示 了 污 
染 物 的 相互 依存 关系 ， 还 揭示 了 哪 一 类 污染 


es。 IOI8。 


源 对 纳 污水 体 的 影响 最 大 ， 即 间接 评价 了 各 
类 污染 源 对 环境 的 污染 程度 ， 将 水 质 与 污染 
源 的 关系 直观 地 表征 了 出 来 。 

3.4 ”从 各 污染 因子 的 相关 和 矩阵 看 出 ,各 污染 
因子 之 间 的 相互 关系 的 网 络 结构 和 层次 结 
构 , 这 对 于 寻找 治理 污染 的 途径 是 很 有 用 的 。 
如 降低 悬 译 物 的 输入 ， 有 可 能 降低 氟 化 物 及 
NOs 一 N、NOs 一 NN 的 输入 ， 从 而 降低 了 
BOD;、COD， 增加 了 DO, 减少 了 NHs 一 N 
的 污染 。 

3.5 从 污染 因子 的 相关 角度 看 , 武 水 宜 章 江 
段 SS 与 NO, 一 N、NO; 一 N、F 、Pb、 Zn 等 
的 相关 系数 在 0. 98 以 上 ,也 许可 表明 水 中 的 
NO 一 N、 NO; 一 N、F 、Pb、 Zn 等 污染 物 的 
来 源 与 面 污染 源 有 密切 关系 。 该 区 有 部 分 区 
域 是 高 所 区 ， 区 域内 有 较 多 小 型 黑色 金属 与 
有 色 金 属 矿 点 ， 这 种 工业 特点 也 进一步 表明 
了 这 一 点 。 


参考 文 献 ( 略 ) 


收 稿 日 期 ，2000 一 05 一 ]2 


